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Abstract: According to recent theoretical developments, it might be possible to 
produce mini black holes in the high energy experiments in the LHC at CERN. 
We propose here a model based on the n-dimensional Vaidya metric in double null 
coordinates for these decaying black holes. The associated quasinormal modes are 
considered. It is shown that only in the very last instants of the evaporation process 
the stationary regime for the quasinormal modes is broken, implying specific power 
spectra for the perturbations around these mini black-holes. From scattered fields 
one could recover, in principle, the black hole parameters as well as the number of 
extra dimensions. The still mysterious final fate of such objects should not alter 
significantly our main conclusions. 
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1. Introduction 



New models with extra dimensions [1] predict the production of mini black holes in 
particle accelerators with sufficient large energies. Such events are expected to be 
obtained in the LHC at CERN[2]. These n-dimensional mini black holes are ex- 
pected to be highly interacting, and, once formed, Hawking radiation [3] is expected 
to settle after possible transient stages. Their phenomenological and observational 
consequences have been intensively discussed. (See [4] for recent reviews). In partic- 
ular, it is expected that their net radiated power, and hence, their mass decreasing 
rate, be driven by the n-dimeusional Stefan-Boltzmann law[5], leading to 
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where tp and Mp stand for the Planck time and mass, respectively, while a„ is the 
effective n-dimensional Stefan-Boltzmann constant [5], which depends on the avail- 
able emission channels for the Hawking radiation[6]. Typically, however, one should 
expect a„ ~ 10~^. Eq. (1.1) can be immediately integrated, 

n — 3 

1 - - ) , (1-2) 
< t < to, where the lifetime to of a black hole with initial mass mo is given by 

n-3 f mo\"--' tp_ 



n-1 \Mp 

Following Arkani-Hamed et al[l], the phenomenology of such mini LHC black holes 
can be studied by setting the Planck scale in order to have Mp ?a ITeV. 

We recall that (1.1) is not expected to be valid in the very final stages of the 
black hole evaporation, where the appearance of new emission channels for Hawking 
radiation can induce changes[6] in the value of the constant a„. Perhaps even the 
usual adiabatic derivation of Hawking radiations is not valid any longer [7]. We do 
not address these points here. We assume that the black hole evaporates obeying 
(1.2) for < i < ^0- The numerical analysis, however, requires the introduction of a 
regularization for the final instants of the evaporation process. Nevertheless, as we 
will see, our main results do not depend on such final details. 

Here, we consider the quasinormal modes (QNM) associated to a radiating n- 
dimensional black hole with the decaying mass (1.2). Since the preferred emission 
channels for Hawking radiation correspond to massless fields, we model these evapo- 
rating mini black holes by means of an n-dimensional Vaidya metric [8, 9] in double- 
null coordinates[10, 11, 12]. The Vaidya metric corresponds to the solution of Ein- 
stein's equations with spherical symmetry in the presence of a radial flow of unpo- 
larized radiation. Such evaporating mini black holes, however, are not expected to 
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emit isotropically[13, 14] and, hence, any realistic model should not be spherically 
symmetric. Furthermore, a typical mini black hole created in the LHC environ- 
ment should not have zero angular momentum. Our simple model, nevertheless, is a 
step toward the construction of more realistic ones. The identification of stationary 
regimes [15] for the QNM in the non-spherically symmetrical case, for instance, could 
simplify the analysis of more realistic configurations. We will return to these issues 
in the last section. 

In n-dimensional double-null coordinates {u,v,9i, . . . , the Vaidya metric 

has the form [9] 

ds'^ = -2g{u, v)dudv + r^{u, v)d^l^_2, (1.4) 

where dVL^_2 stands for the unity (n — 2) dimensional sphere, and g{u, v) and r(M, v) 
are smooth non vanishing functions obeying[12] 

g = -9„r, (1.5) 

^^'^2- (n-3)r--3 ' ^^"^^ 

where, for the present case of an outgoing radiation field {w! < 0), corresponds to 
the retarded time coordinate. We adopt hereafter natural unities (tp = Mp — — 

1)- 

Our choice for the mass function m{y) is guided by the solution (1.2). Never- 
theless, the final stage of a black hole evaporation is a rather subtle issue. A black 
hole could evaporate up to zero mass as described by (1.2), leaving behind an empty 
Minkowski-like spacetime[16] (or even a naked singularity[17]), or it could evaporate 
until it reaches a minimum mass, that is, turning into a massive remnant [18]. In order 
to circumvent these problems in the numerical analysis, we introduce a regularization 
for the final stage of the evaporation process. We consider the mass function 

{ n-3 
A — Stanhp(i> — i>i), V > 

with p > 0. The constants A, B and p are determined by imposing conditions for 
the continuity of m[v) and its first derivative in v = v\. Clearly, A — B = mp 
is the mass of the final remnant. The regularization is efi^ective only at the final 
instants of the evaporation process, (fg — v^jv^ ^ 1, and mp <^ ITeV. We will 
show that, during the major part of the evaporation process, the stationary regime 
for the quasinormal modes described in [15] holds, implying a specific, and perhaps 
observable, power spectrum for the perturbations around these mini black-holes. 
Besides, the perturbation power spectra should not depend significantly on the final 
fate of the black-hole. 
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2. The Quasinormal modes 

We decompose a generic perturbation field 4> as 

^ = J]r-^V'^(ii, v)Yem{ei, 6^-2), (2.1) 

im 

where stands for tlie scalar spherical harmonics on the {n — 2) unity sphere, for 
which <9^^_2^m ~ —i{£ + n — 3)1^^, where £ = 0, 1, 2, . . ., and m denotes a set of 
(n — 3) integers (mi, 1712, . . . , rrin-s) satisfying i > nins > m2 > \mi\ (See [19] for a 
concise description of higher dimensional spherical harmonics). By using (1.5) and 
(1.6), the equation for ip^ reads 

^ + giu,v)Viu,v)i.e-0, (2.2) 

where 

.(„..)4(!(i±pl),(!i^||i^,(i^^^ (2.3) 

The constant a determines the type of the perturbation considered: a — cor- 
responds to scalar and gravitational tensor perturbations, cr = 2 to gravitational 
vector perturbations, a — 2/{n — 2) to electromagnetic vector perturbations, and 
(7 = 2 — 2/(n — 2), finally, to electromagnetic scalar perturbations[20]. 

We perform an exhaustive numerical analysis of the equations (1.5), (1.6), and 
(2.2) along the same lines of the method proposed in [15]. In particular, we could 
verify that the QNM stationary behavior for slowly varying masses reported in [15] 
is not altered in higher dimensional spacetimes, see Fig. 1. Hence, provided the 
mass function m{v) varies slowly, the QNM of (2.2) set down in a stationary regime, 
and the associated frequencies (cDr) and damping terms {uji) follow the track corre- 
sponding to l/rh{v), where Vh is the (aparent) horizon[ll] radius of an n-dimensional 
black- hole of mass m{v). In a more quantitative way, one has here 

C^R,l(^) _ ^fe(O) _ A _ _^ 

where a;Rj stand for the oscillation frequency (r) and damping term (i) of the QNM 
corresponding to an n-dimensional Schwarzschild black hole with mass mo = m(0). 
We notice that the relation c<jr cc l/r^ for n-dimensional Schwarzschild black holes 
has been previously obtained by Konoplya in [21]. As in [15], we have used Gaussian 
initial conditions for all the analysis, although equivalent results can be obtained for 
any localized initial condition. 

Our simulations strongly indicate that the condition for the QNM stationary 
regime[15] must be generahzed for the case of n-dimensional black holes as |r^'(T;)| < 
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Figure 1: QNM "instantaneous" frequencies (see [15]) for the equation (2.2). In the 
stationary regime, the frequencies follow the track corresponding to l/r}i{v). The cases 
depicted here corresponds to a = 0, £ = 2, and a = 0.02. Note that, the smaller the value 
of Qn, the better the fit. 



|ci;i(t;)|, where uji{v) is the smallest damping term of the system. In the present case, 
it reads 

21 



1-^ 



1 71-2 

'n-1 



n 



n 



(2.5) 



where uoi is the smallest damping term of an n-dimensional Schwarzschild black 
hole with mass mo = m(0), typically corresponding to scalar perturbations. For 
the LHC mini black holes, the term between square brackets should be of order 
of unity, irrespective of n. Thus, only in the very late times of the evaporation 
process (for typical small values of a„, for less than the last ai""^-*^^""^^ fraction of the 
lifetime period) the stationary regime is broken. Hence, the late time exponentially 
suppressed perturbation of (2.2) is well approximated by 



^.(t,) = e-'^'^ sin ((Dr^; + (5) 



(2.6) 



for < < ?;0) and ^(f ) = for f > fo, where (Drj are themselves functions of v 
given by (2.4), and 5 is an arbitrary phase. 

Our main observation is that, for the typical values of the parameter a„ and 
mo ~ ITeV, the Fourier spectrum ^(/) of the stationary perturbations (2.6) is very 
close to the the Fourier spectrum ^(/) of the perturbations corresponding to a n- 
dimensional Schwarzchild case {m{v) — mo). 



'0(v) = e "^'^ sin {uj-rv + 5) 



(2.7) 



for V > This fact, clearly illustrated in Fig. 2, certainly would deserve a more 
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Figure 2: The power spectra: |^(/)|, the sohd hne, given by Eq. (2.12); and |^'(/)|, 
the hne with points, calculated numerically from (2.6) and (2.4). They arc indeed very 
close for the typical values of a and tuq ~ ITeV. In particular, both spectra exhibit similar 
pronounced peaks. Note that the discrepancies for large values of u; are due to the FFT 
aliasing effect for frequencies larger than the Nyquist critical frequency [22], and not to 
real differences between |^'(/)| and |*(/)|. The case depicted here corresponds to n = 4, 
04 = 0.002, ujR = 0.25, ui = 0.01, and 5 = 0. 



rigorous analysis. Some simple estimations, however, do endorse the observation. 
From the linearity of the Fourier transform and Parseval's theorem, we have that 



dv — 



^if)-^{f) df. 



{2.i 



Hence, provided the left handed side of (2.8) be small, ^{f) will be close (in the 
norm) to ^(/). The integral in the left handed side of (2.8) can be split as 



h + h 



dv + 



V2 



^-0 — ijj^ dv. 



The second integral can be estimated as 



h<2 



( I iP'^dv + I i^^dvj < 4 /" 



-2u}\V2 



e-'^^^'dv = 2- 

V2 



(2.9) 



(2.10) 



Typically, a;R, and LiJ\ are of the same order (unity), while a„ is much smaller (10~^). 
If one chooses V2 corresponding, for instance, to 10 oscillation cycles of '?/'(f), the 
value of I2 will be less than e^"^^ . This is the error involved in approximating the left 
handed side of (2.8) by l\. On the other hand, during the 10 first oscillation cycles 
of '?/'(f), the variations of d;R^i(f) are of the order 1 — (100/99)^/'-"^^-' for black holes 
with initial mass mo = 1 TeV. Hence, in the interval [0, ^2], t/'(w) is indeed very close 
to ■0(f) r implying that the left handed side of (2.8) is small and, finally, that ^(/) is 
close to *(/). 
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Figure 3: Graphical solution of (2.13). It is clear that, for arbitrary phases 5, the peak 
of the power spectrum (2.12) is located in the range given by (2.14). 



The Fourier spectrum ^(/) of the perturbation (2.7) can be easily calculated 

^Jo^' ^ (a;R)' + (^/ + WI)' ^ ^ 

The associated power spectrum 



. (cjR COS ^ + cji sin + P sin^ ^ 
has a pronounced peak (see Fig. 2) at jmax given by 



from where we can conclude that 



-uJl<j^^<^pi^l (2.14) 
provided |a;R| > |a;i|, see Fig. 3. 



3. Final Remarks 

The most interesting part of our results concerns the characterization of the signals 
that might come out of the black hole probe. From the peaks in the power spectra of 
the perturbations around the evaporating mini black holes, it is possible to determine 
uj-^ and LO\ and, consequently, infer some of the black hole parameters as its initial 
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mass mo and even the dimension n of the spacetime where it effectively hves. We do 
not expect, of course, the gravitational perturbations associated to these mini black 
hole to be measurable. However, we remind that the QNM analysis can be applied for 
any test field propagating around the black hole. In particular, it also applies for real 
(confined to the brane) electromagnetic perturbations. Moreover, since these higher- 
dimensional mini-black holes are expected to emit large fractions of their masses 
in the bulk [13], the assumption of a spherically symmetrical n-dimensional Vaidya 
metric can be justified first approximation. 

According to the discussion of last section, the late time behavior of electromag- 
netic waves scattered by these evaporating mini black holes should exhibit a power 
spectrum as that one depicted in Fig. 2, since the electromagnetic perturbations will 
be also of the form (2.6) for large times. Furthermore, the energy carried by these 
perturbations is negligible and, thus, they are expected to be fully sensitive to the 
higher dimensional dynamics of the black hole. Hence, from a precise determination 
of the peak location for real electromagnetic perturbations, we can get the relevant 
parameters of the mini black hole, including the number of extra dimensions. For 
instance, for 4-dimensional black holes, the QNM frequencies and damping terms 
for the first electromagnetic perturbations [23] are cjr = 0.2483 and Ui = 0.0925, 
implying that the frequency peak of the electromagnetic power spectrum be in the 
range 



Typically, the larger is the number of extra dimensions, the larger will be the peak 
frequency, even surpassing the limit of 1 TeV. However, one can, in principle, de- 
termine the relevant parameters of the mini black hole by comparing the peaks 
of power spectra of real (confined to the brane) perturbation fields with the calcu- 
lated higher-dimensional QNM frequencies and damping terms for the n-dimensional 
Schwarzschild black-hole [24]. 
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